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ABSTRACT 



An experimental Investigation of the internal ballistics 
of solid fuel ramjets was conducted in order to determine 
the regression rate of the fuel as a function of chamber 
pressure, inlet air temperature, and air flux rate, and to 
model the flow in solid fuel ramjets which use sudden expan- 
sion flameholders at the inlet. In addition, flame stabili- 
\ 

zation limits were investigated. A computer solution for 
the non-reacting flow field gave results in good agreement 
with experiments. 

Solid fuel ramjets have an average regression rate of 
the fuel that closely follows the theoretical expression 
derived for klnetically controlled hybrid rocket combustion. 
The inlet step-height-to-motor-diameter ratio is the domi- 
nant parameter in determining flame stability limits. 
However, inlet velocity is also an important parameter. 
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I. INTRODUCTION 



The solid fuel ramjet can be classified as a special 
case of the hybrid rocket. That is, the fuel is solid and 
the oxidizer is gaseous (air). The solid fuel ramjet pro- 
vides many of the advantages of solid fuel rockets without 
the weight penalty of carrying its ovm oxidizer. Ease of 
handling, storability, weight savings, and possibly greatly 
lncrease(i tactical ranges are among its advantages. 

A research project is currently being conducted at the 
Naval Postgraduate School to model the internal ballistics 
of solid fuel ramjets. Such a model would provide the 
design engineer with a working tool that would permit 
performance parameters and operating characteristics to be 
predicted. The overall goal of the research project is to 
develop a model that can predict the regression rate of the 
fuel as a function of operating and configuration variables. 
Very little experimental data currently exist that can be 
used to determine the adequacy of a model, thus the model 
must be developed in steps verifying each step with experi- 
mental data. Non-burning tests to determine flow patterns, 
non-burning tests with wall blowing to determine the effect 
of wall mass addition on the flow characteristics, and hot 
firing tests with various inlet geometries are required. 

This thesis initiates the overall research project 
obj ectlves . 
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Information is provided from several areas of study 
that are pertinent to the internal ballistics of solid fuel 
ramjets. Specific areas are flov/ and heat transfer in. ducts 
with a backward facing step, flame stabilization in high 
speed flows, hybrid rocket combustion, and thermal decompo- 
sition of polymers. 

Abbott and Kline [Ref. 1] investigated the turbulent 

flow over backward facing steps using water flow through a 

\ 

two-dimensional channel. It v;as found that the recircula- 
tion zone consisted of a complex pattern involving three 
distinct regions. Immediately downstream of the step face 
there is a three-dimensional zone "characterized by one or 
more vortexes rotating about an axis normal to the vertical 
wall." This zone is followed by an almost two-dimensional 
zone of recirculating flow. The area at the flow reattach- 
ment point is characterized by a fluctuating three- 
dimensional flow. It was also found that the location of 
the reattachment point remained fixed as Reynolds number 
and/or inlet turbulence intensity varied, and that increas- 
ing the step height moved the reattachment point further 
downstream from the step face. 

Krall and Sparrow [Ref. 2] using water flow through an 
orifice in an electrically heated circular tube found that 
local heat transfer coefficients in the separated, reattached, 
and redevelopment regions are several times larger than those 
for fully developed turbulent pipe flow. They found that 
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the maximum heat transfer coefficient occurred in the 
vicinity of the point of flow reattachment. 

At low Reynolds numbers based on downstream pipe . 
diameter (10,000 to 50,000) the heat transfer near the 
reattachment point varied rapidly with Reynolds number. 
However, at higher Reynolds numbers (70,000 to 130,000) 
only a weak dependence was noted. They also found that for 

weaker separations (small step heights, h) the reattachment 

\ 

point was' well defined but for strong separations (large 
step heights) the reattachment point was spread out. This 
transition from a single well-defined reattachment point to 
the "spread-out" configuration occurred between h/D=0.l6 
and h/D=0.25, where D is the downstream pipe diameter. The 
reattachment point was located between 1.25 and 2.5 pipe 
diameters downstream of the Initial separation. Increased 
step height moved the reattachment point slightly downstream. 
Reynolds number and Prandtl number had no significant effect 
on the reattachment point location. 

The results of Krall and Sparrow are also in general 
agreement with the results of Abbott and Kline. For small 
step heights (h/D=.17) the two-dimensional and axisymmetrlc 
reattachment points are very nearly equal (X^/D from 1.0 to 
1.25). However, for larger step heights, the reattachment 
zone spreads out considerably more for the two-dimensional 
case than for the axisymmetrlc case. 

The ideal hybrid rocket is one in which a gaseous 
oxidizer flows over a solid fuel. In a hybrid rocket, gas 
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phase combustion in the form of a diffusion flame takes 
place in the turbulent boundary layer formed by the oxidizer 
flow over the solid surface of the fuel. Heat is trans- 
ported by convection and radiation to the solid surface 
which causes decomposition of the fuel. The convective heat 
transport is complicated by the presence of high rates of 
wall mass addition. Using a heat transfer limited model; 

regression rates of the fuel are generally proportional to 

\ 

G^, where G is the total gas flux and n is a constant 
(typically 0.8). 

Experimental data obtained at low chamber pressure 
and/or high oxidizer mass flux indicate that the regression 
rate of hybrid grains becomes less dependent on G, but is 
directly related to the partial pressure of the oxidizer in 
the gas phase [Refs. 3 and Reduced reaction rates at 

lower pressures broaden the flame zone. Schlieren studies 
have also shown that the flame zone broadens at low chamber 
pressures and may broaden enough to undergo transition to 
a turbulent premixed flame [Ref. 6]. When the oxidizer is 
diluted, the diffusion flame moves further away from the 
surface, increases in thickness, and the flame temperature 
drops [Refs. 3 and 53. In addition, surface kinetics may 
be rate controlling under these operating conditions 
[Ref. 7]. When these conditions exist, the heat transfer 
limited model does not adequately describe the combustion 
process. 
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Hybrid rockets operating at high chamber pressures are 
adequately described by a heat transfer limited model. 

Muzzy [Ref. 4] has revlev;ed the modifications that have been 
made to the heat transfer limited model to account for gas 
phase and surface kinetic effects. The modified models 
[Refs. 6 and 7] predict a dependence of the regression rate 
on pressure as well as mass flow rate. For a system where 

radiative heat transfer can be neglected. Muzzy has shown 

\ 

that the regression rate equation In the klnetlcally domi- 
nated region of operation can be expressed as: 

r = C 

where C Is a constant, P Is pressure, G Is total mass flux, 

X Is the longitudinal coordinate, a Is a constant that 
varies from 0.3 to 0.^, and b Is a constant that varies from 
-0.1 to -0.2. Another study by Kustov and Rybanln [Ref. l8] 
has also treated the pressure sensitive regime of operation. 

In a solid fuel ramjet the reacting flow has character- 
istics that are similar to those discussed above. However, 
two significant differences exist. Solid fuel ramjets would 
normally operate at lower chamber pressures than hybrid 

r 

rockets, resulting In significantly greater port velocities. 
In addition, the oxidizer In this case is air which would 
correspond to a severely diluted oxidizer in hybrid combus- 
tion. These differences require some means of flame stabi- 
lization. A series of tests performed at the United 
Technology Center [Ref. 8] have shown that one means of 
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sustaining combustion in a solid fuel ramjet is to use a 
sudden expansion inlet as is used in dump combustors. 

Similar results were reported by McCarthy [Ref. 9]. 

In a high velocity gas stream, where the velocities of 
the flow are at least an order of magnitude greater than the 
turbulent flame speeds of the combustion mixtures, flame 
stabilization is required. Bluff body flame stabilization 

has been used successfully in turbojet and liquid fuel 

\ 

ramjet combustor design in the past. Typical flame stabi- 
lizers have been V-gutters and rods, or the dome of a can- 
type combustor [Refs. 10 and 11]. In addition, dump or 
sudden expansion combustors have been studied for Integral 
rocket/liquid ramjet designs [Ref. 19]* 

The recirculation zone behind the bluff body or step 
causes Intense mixing of fresh unburned fuel and air with 
hot combustion products. The hot combustion products ignite 
the fuel-air mixture. This stabilized flame can then 
propagate throughout the combustor. Due to the intense 
turbulent mixing in the recirculation zone, combustion may 
in some cases approach that of a well-stirred reactor. 

A well-stirred reactor model considers the composition 
and all of the thermodynamic properties to be uniform 
throughout the volume. An idealized reactor with instanta- 
neous mixing would have homogeneous fuel, air, and product 
concentrations, and constant pressure and temperature 
[Refs. 12 and 13]. 
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Solid fuel ramjet combustors have, depending upon the 
combustor length and inlet configuration, two separate and 
distinct combustion zones, each of which appears to be 
controlled by a different mechanism [Ref. 9] (See Figure 1). 
In the recirculation zone products and reactants are rapidly 
mixed, and in the limiting case a well-stirred reactor may 
be approached. Downstream of the flow reattachment point 

the boundary layer develops and the combustion may be similar 

\ 

to that in a hybrid rocket in which kinetic effects are 
dominant (low pressure, high mass flux, severe oxidizer 
dilution) . 

Thus, one possible model for solid fuel ramjet combus- 
tion would be to treat the recirculation region as a well- 
stirred reactor and the downstream region as a kinetically 
dominated hybrid. The larger the grain length-to-dlameter 
ratio, the more dominant the downstream region becomes in 
determining the average regression rate, and the more one 
could expect the ramjet to behave as a kinetically controlled 
hybrid combustor. However, for small L/D ratios and/or 
large recirculation zones, the combustion in the recircula- 
tion region may become important in determination of the 
average regression rate. In addition, the combustion char- 
acteristics in the recirculation region determine the 
flammability limits for the entire combustor. Thus, the 
model used in this region requires considerable attention. 

For very strong recirculation zones, a well-stirred reactor 
model may be adequate. However, for weaker recirculation 
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zones, which are desirable from a pressure loss standpoint, 
kinetic rates may be of the same order or greater than the 
reactant mixing rates. In this case, a well-stirred reactor 
model cannot be expected to be applicable and a more detailed 
consideration of the recirculating flov; Is necessary. Thus, 
the work of Gosman, et al. [Ref. l 4 ] and Spalding, et al. 
[Ref. 15] becomes Increasingly Important and may provide 
the basis for modeling the recirculation region In solid 
fuel ramjet combustors. 

It Is the purpose of this study to determine the 
dependence of the regression rate of a solid fuel ramjet 
on chamber pressure. Inlet air temperature, and mass flux 
of air, and to model the flow In the recirculation region 
caused by the sudden expansion Inlet. 
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II. METHOD OF INVESTIGATION 



Non-burning flow visualization tests vrere conducted 
using a plexiglass model with three inlet step configura- 
tions. Flov; reattachment points of the incoming stream 
were measured for various step Reynolds number flows. 

A Fortran language computer program, Plstep II, devel- 
oped by Spalding, et al. [Refs. and 15] provided a finite 
difference solution of the two dimensional Navier-Stokes 
equations governing heat and mass transfer in recirculating 
flows. It was used to calculate the flow reattachment point 
and the flow field for non-burning flows. The results were 
compared with the flow visualization data and the work of 
previous investigators. The boundary conditions on the 
differential equations were altered to simulate' blowing at 
the wall in order to study the effect of mass injection 
from the wall on the flov; field. 

Experimental firings were conducted using polymethyl- 
methacrylate grains burning with non-vitiated air to provide 

regression rate data as a function of inlet air temperature, 

( 

chamber pressure, and oxidizer flux rate. These tests were 
made to provide initial data for characterizing ramjet 
combustion and flame blow-off limits, and to determine if 
the reacting flow behavior was similar to that obtained 
analytically and experimentally for non-burning flow. 
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III. DESCRIPTION OP APPARATUS 



A. COLD FLOW VISUALIZATION MODEL 

The flov; visualization model (Figure 2) consisted of 
three sections. The first section was a 2% x 2^5 x 5-inch 
plexiglass block with a three-quarter inch diameter hole 
bored through the center. It received air from high pressure 
tanks and supported a perforated disk flow straightener . 

The step section was made of stainless steel and was ten 
inches in length to further dampen turbulence and provide 
fully developed velocity profiles. Three steps were used 
with inside diameters of 0.246, 0.519, and 0.714 Inches 
respectfully. The third section was a 15%-inch plexiglass 
tube with an inside diameter of 1.469 inches. 

For the lower flow rates a sonic choke with a diameter 
of 0.08l Inches was used to measure flow rates. For higher 
flow rates, an orifice flow meter with inputs to a Vlsicorder 
was used to measure flow rates. The flow reattachment 
points were measured visually using a scale. 

B. HOT FIRING APPARATUS 

A schematic of the solid fuel ramjet apparatus showing 
the overall system is presented in Figure 3- 

1 . Ramjet Motor 

The solid fuel ramjet motor consisted of four main 
sections; the head end assembly, the step insert section, 
the grain, and the aft closure section (Figure 4). All of 



17 



the metal parts of the motor except for the bolts were made 
of stainless steel. 

The head end assembly was 6h Inches long and contained 
the inlets for the air, the oxygen-methane ignition mixture, 
the nitrogen purge, and the cooling air. Two Champion Z7 
spark plugs v:ere used for ignition. 

The step insert section was 2-3/4 inches long and 
held the step in place. Three steps were used during this 
investigation with inside diameters of 0.25, 0.50, and 0.75 
inches. This section also accomodated the forward end of 
the PMM grain. 

Polymethylmethacrylate was chosen as the fuel in 
this study because it has been widely used for basic research 
in hybrid rocket combustion and in studies of polymer 
degradation in the past. It is composed of long chains of 
a repeating monomer: 

CH^ 

C 

COOCH^ 

Idealized chemical reaction with complete combustion 
would involve a stoichiometric reaction with air as follows: 

^5^8°2 22.56N2 ^ 4H2O + 5CO2 + 22.56N2 

The average molecular weight of the products in this case 




18 



is 31*3 Ibm/lbmole. The gas constant is R/M = 
ft-lbf/lbm-°R, and the stoichiometric air fuel ratio is 
8.27 Ibmair/lbmfuel . The equilibrium ratio of specific heats 
at a temperature of 3500°R is y = 1.24. An adiabatic flame 
temperature calculation for the stoichiometric reaction using 
a value of of 92.8 kcal/mole for the heat of formation of 
PMM at 25 °C [Ref. 5] gives a flame temperature (with no 
dissociation) of 4265 °R. 

The polymethylmethacrylate grains were cut from 
a slab into 3-3/4 x 3-3/4 x 12 inch blocks. The ends were 
rounded to an outside diameter of 3-1/2 inches to fit into 
the step section and the aft closure section, and were 
sealed with "0" ring seals. The Initial inside port diameter 
of the grains was 1.5 inches. Eight grains were made that 
were six inches in length, but the majority of the tests 
were performed using grains that were 12 Inches in length. 

The aft closure section contained a straight 
section with a pressure tap and a converging-diverging 
nozzle. Nine nozzles with various throat areas were used 
to control the chamber pressure. A Wiancko 0-200 pslg pres- \ 

V 

sure transducer was mounted on the forward part of the aft 
closure section to measure the chamber pressure. 

The grains were held in place by four 3/8 inch 
diameter steel rods between the step section and the aft 
closure section. Mounting flanges on the head end assembly 
and the aft closure section were bolted to a large test 
stand. 
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2 . Ignition System 



The Ignition system metered oxygen and methane at 
800 pslg from two high pressure tanks through needle valves 
Into the head end assembly section. Inlet ports on the top 
and the bottom Injected methane, and two Inlet ports 
injected oxygen from either side. Two Inches downstream 
from the oxygen-methane ports were two Champion Z7 spark 
plugs located l80 degrees apart. A transformer supplied 
power to a Model T spark coll to energize the spark plugs. 

3. Nitrogen Purge and Cooling Air System 

The nitrogen purge and cooling air system was 
connected upstream of the head assembly section. Two high 
pressure tanks fed nitrogen into the motor at the end of a 
hot firing for approximately one second to extinguish the 
combustion flame. Cooling air from a low pressure compressor 
was then blown through the motor. 

4 . Air Flovf Control 

A standard ASME orifice flow meter was used to 
measure the flow rate of the air into the motor. A manually 
operated gate valve between the orifice and the motor was 
used to provide the desired flow rate through the motor. 

Two pneumatically operated Jamesbury ball valves operating 
together either vented the air to the atmosphere, or allowed 
the air to pass through the motor. 

5. Data Acquisition System 

A Colvin 0-35 psl differential pressure transducer 
was used to measure the pressure drop across the orifice 
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flow meter. A Wlancko 0-200 psig pressure transducer was 
used to measure the pressure upstream of the orifice, and 
an Iron-constantan thermocouple v;as used to measure the 
temperature at the orifice. A second iron-constantan thermo- 
couple was used to measure the temperature of the air enter- 
ing the motor. A Wiancko 0-200 psig pressure transducer 
measured chamber pressure at the aft closure section of the 
motor. 

\ 

The orifice thermocouple output was recorded on a 
0-600°F strip chart recorder. All of the other transducer 
outputs were recorded on a Honeywell Model 2106 Visicorder. 

A five cycle per second timing signal and an ignition pulse 
were also recorded on the Visicorder. 

6 . Air Feed System 

A Pennsylvania air compressor provided air at pres- 
sures up to 150 psia at a rate of up to 700 scfm. This air 
was fed into an air reservoir and a Polytherm air heater. 

The Polytherm air heater was capable of yielding 1.75 pounds 
per second airflow at 150 psia and 1000°F. The heater 
burned natural gas, and the air passed through a heat 
exchanger, thus providing non-vitiated hot air. 

A temperature controller directed two mixing valves 
that produced the required air temperature at 150 psia. A 
manual gate valve controlled the air flow to the test cell. 

A valve in the test cell controlled the flow of air to the 
pipe leading to the orifice flow meter. 

Figures 5> 6, 7, and 8 show the test cell and the 
ramjet apparatus. 
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IV. EXPERIMENTAL PROCEDURE 



A. COLD FLOW TESTS 

The cold flow tests were performed at the rocket lab 
facilities of the Naval Postgraduate School. Visualization 
was obtained by injecting a mixture of water, alcohol, and 
food coloring at the base of the step inlet through a small 

tube with a hypodermic needle. The mixture was broken up 

\ 

into very small droplets by the turbulence in the recircu- 
lation zone. In the recirculation zone the droplets on the 
wall moved towards the step. At the flow reattachment point 
the droplets on the wall were stationary. Downstream of the 
reattachment point the droplets on the wall moved away from 
the step. Using this technique, the flow reattachement 
point locations for various step Reynolds number flows using 
three step sizes were measured visually using a scale. 

Initial visualization tests were attempted using a 
tufted model but they proved to be of limited value due to 
the highly turbulent nature of the flow. 

B. HOT FIRING TESTS 

All hot firing tests were performed in the jet engine 
test cell at the Naval Postgraduate School. Hot firing 
data were obtained from a series of firings at three air 
flow rates, three chamber pressures, and three inlet temper- 
atures, using three step heights and two grain lengths. The 
majority of the hot firings were performed with a 12 inch 
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grain length and an inlet diameter of 0.5 inches. Table I 
summarizes the nominal test conditions. 



TABLE I 

NOMINAL TEST CONDITIONS 



Chamber Pressure (psia) 


50 


75 


100 


Air ^low Rate (Ibm/sec) 


.15 


.25 


.40 


Inlet Temperature (®P) 


70 


200 


350 


Inlet Diameter (in) 


.25 


.50 


.75 


Step h/D 


.i|17 


.333 


.25 


Grain Length (in) 


12 


6 





For the tests in which the Polytherm air heater was used, 
the temperature controller was set and the temperature would 
stabilize at the ramjet inlet in one to two hours. This 
delay was due to the fact that the flow rates to the motor 
were quite low, at most 0.5 Ibm/sec. Once the temperature 
stabilized at the desired temperature at the inlet to the 
ramjet, the manual gate valve downstream of the orifice was 
set to obtain the desired air flow rate. 

Chamber pressures were preselected and controlled by 
varying the nozzle throat area for a given air flow rate. 
Chamber pressures could normally be attained within ten psi 
of the desired pressure. Calculated combustion temperatures 
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varied significantly with air flow rate and step size. Also, 
at higher chamber pressures the preset flow rate would 
decrease upon ignition because the higher back pressure would 
unchoke the air flow control valve. 

Ignition normally lasted for four seconds. After three 
seconds of ignition, the air was directed through the ramjet 
motor. Ignition was continued into the first second of the 

run to Insure that combustion would be sustained. A series 

\ 

of thirteen tests at high flow rates in which the ignition 
flame was blown out \^ere used to determine the rate of 
consumption of the PMM grains during the oxygen-methane 
Ignition. These data were used to correct the initial 
weight of the fuel used in the regression rate calculation. 

Combustion lasted for thirty seconds. At the end of 
each run the motor was extinguished by simultaneously venting 
the air to the atmosphere with the Jamesbury ball valves and 
actuating the nitrogen purge system. Low pressure air v;as 
then blown through the motor for cooling. The only section 
of the motor that heated up significantly was the aft 
closure section and the nozzle. 



V. PISTEP II COMPUTER PROGRAM 



Pistep II is a Fortran computer program developed by 
Gosman, Pun, Runchal, Spalding, and Wolfshtein [Ref. l4]. 

It provides a finite difference solution of a steady 
incompressible turbulent flow through a cylindrical pipe 
with a sudden step enlargement. Equations [Ref. 15] 

governing the distribution of stream function, ip, vorticity, 

\ 

co/r, non-dimensional temperature, T, turbulent kinetic energy, 
k, and turbulent energy dissipation rate, e, are cast in 
elliptic form in cylindrical polar coordinates as follows: 




9z 



If) - I? If) ^ 



9z 






3(C^4.) 

9z 



] 








3r 






0 



where <f) represents any one of the variables w/r, k, e, 
or T. The functions a^, b^, C^, and are presented in 
Table II. 

A two equation turbulence model is employed in the 
program. The two variables are turbulent kinetic energy, k, 
and the turbulent energy dissipation rate, e. The effective 
viscosity, is related to k and e by 



^eff 



C pk^/e 

y 
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A "dissipation length scale", l, can be used in an 



alternative formulation: 



“eff = 



Equating the two equations for the effective viscosity gives 






Constants that appear in the two turbulence equations 
are c^, 02 » c^, Cj^, and a^. Values used in the computer 
program together with their Fortran symbols are given in 
Table III. 



TABLE II 

VARIABLE FUNCTIONS (Table I of Ref. 15) 



4- 








^4* 


o)/r 


r2 


2 ' 
r 


^eff 


0 


’I' 


0 


1/pr^ 


1 


-o)/r 


k 


1 


^eff^^k 


1 


-(G-pe ) 


£ 


1 


^eff^^e 


1 


2 

-(c^Ge/k - C2pe /k) 


T 


1 


^eff/% 


1 


0 


G is 
G = 


a generation term: 

2 ^^r 


■)^ + 


V^ o 3V 8V^ p 



26 






I 



TABLE III 



TURBULENCE CONSTANTS (Table II of Ref. 15) 



Constant 


Cn 


Co 


c 




0 


1 


2 


y 


k 


e 


Fortran 

Symbol 


Cl 


C2 


CMU 


PR(NK) 


PR(NE) 


Value 


1.45 


2.0 


0.09 


1.0 


1.3 


Boundary 


conditions 


are 


specified 


at every 


point on the 



boundary. Wall functions that are derived from Couette flow 
solutions are used in evaluating boundary values for w/r, k, 
e, and T. Other boundary conditions are given in Figure 9. 

The computer program calculates heat and momentum 
transfer in steady two-dimensional axisymmetric recirculating 
flows. Five simultaneous coupled non-linear partial differen- 
tial equations are solved. The finite difference equations 
are Implicit equations that are coupled and non-linear. The 
Gauss-Seidel point iterative method is used to solve the 
finite difference equations. Provisions are made for both 
overrelaxation and underrelaxation of the variables. 

References l4 and 15 contain a more complete discussion 
of the mathematical formulation of the problem, development 
of the finite difference equations, and solution algorithms. 

An excellent discussion of several turbulence models can be 
found in Reference l6. 
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The boundary conditions on the pipe wall can be changed 
so that the effect of mass addition through the wall for 
non-reacting flows can be investigated analytically. A 
constant radial velocity at the wall, which yields a 
linearly varying stream function along the wall, is used to 
simulate wall mass injection (See Figure 10). In this 
initial study, the main flow and the wall mass addition 
were Identical in chemical composition and temperature . 

In Plstep II, wall functions based on Couette flow are 
used to specify the boundary conditions on the pipe wall 
for all ‘of the variables except stream function. For very 
low blowing rates, which are characteristic of solid fuel 
ramjets, Couette flow solutions may yield realistic results. 
However, for higher mass injection rates, Couette flow 
solutions may not adequately represent the boundary 
conditions . 

The centimeter-gram-second system of units is used in 
the Plstep II computer program. 
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VI. RESULTS AND DISCUSSION 



A. COLD FLOW EXPERIMENTS 

A series of cold flow visualization tests were conducted 
to determine the reattachment point location for flow over a 
rearward facing step at the inlet of a cylindrical tube. 

Defining a step Reynolds number as: 

T, p V h ^ 

Re„. = ^ where 

step p 



p... fluid density 

V. . . average velocity in the port of the tube 
h. . . step height 
y... fluid viscosity 



it was found that the reattachment point location remained 
constant as the step Reynolds number was changed for a given 
step height. The reattachment point moved further downstream 
with increasing step height. The tube inside diameter was 
used to non-dimensionallze the reattachment point location 
from the face of the step. A linear least-square curve fit 
of 47 data points was used to obtain the following equation 
(See Figure 11 ). 



X 

= 7.15 (h/D) + 0.35 

Figure 12 compares the data obtained in this investigation 
with the data of Abbott and Kline [Ref. 1 ] and Krall and 
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Sparrow [Ref. 2]. The unsymmetrlc reattachment points 
obtained by Abbott and Kline are shown. The flow reattached 
at different locations on either side of the two-dimensional 
channel. The location of the maximum heat transfer coeffi- 
cient found by Krall and Sparrow for water flow through an 
orifice in a cylindrical tube was used to locate the 
reattachment point data shown in Figure 12. 

The data from this study fell between that of Krall and 

\ 

Sparrow dnd the downstream reattachment point found by 
Abbott and Kline. Air was used in this study whereas water 
was used in the other studies; however, each study found 
that reattachment point location was not affected by changing 
Reynolds number for a constant step height. Krall and 
Sparrow also found that changing Prandtl number of the flow 
did not change the reattachment point location. 

The spread in the cold flow data in Fig. 11 is probably 
due to larger relative errors in measurement closer to the 
step face and not due to Reynolds number effects. 

At higher step Reynolds numbers small vortices were 
observed at the step face similar to the observations of 
Abbott and Kline for two-dimensional channel flow. The 
reattachment point was axially symmetric for most Reynolds 
number flows.. However, at high Reynolds numbers the reattach- 
ment point varied around the circumference of the tube. Four 
distinct lobes could be seen when this occurred and may be 
associated with the small vortices at the step face. It 
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was also observed that the rotation and mixing of the flow 
in the recirculation zone became more Intense as the height 
of the step increased. 

B. PISTEP II RESULTS 

The Plstep II computer program was run for four step 
configurations and five inlet velocities. Table IV summarizes 
the reattachment point results. A linear least-square curve 
fit is shown in Figure 11 and Figure 12 for comparison with 
the cold flow data and the data from other Investigations. 

From these figures it is seen that Pistep II gives a depen- 
dence of X^/D on h/D which is in general agreement with 
experimental results. With the current turbulence model and 
wall functions, the results agree closely with those of 
Krall and Sparrow and fall somewhat below the data obtained 
in this investigation. 

A 21 X 21 grid was used for the h/D of .25 and a grid 
of 20 X 20 was used for the other values of h/D. The 
computing time required for each run was on the order of 20 
to 30 minutes. Underrelaxation of the variables was required 
to ensure numerical stability. Even with underrelaxation, 
the case with h/D of 0.i<21 was numerically unstable for 
most of the inlet velocities tried. However, this case 
converged the fastest when numerically stable. 

At each of the grid points, Pistep II calculated the 
distribution of the following parameters: vorticity, stream 

function, turbulent kinetic energy, turbulence dissipation 
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TABLE IV 



RESULTS OF PISTEP II CALCULATIONS 



h/D 


Xj,/D 


Re . * 

step 


Niter** 


VfnC cm/sec) 


.250 


1.735 


2446 


900 


1524 


.250 


1.696 


12232 


700 


7620 


.250 


1.761 


24465 


900 


15240 


.250 


1.702 


36698 


700 


22860 


.289 


1.979 


2009 


800 


1524 


.289 


2.008 


20089 


800 


15240 


.289 


2.004 


31033 


800 


22860 


.3^2 


2.152 


1335 


700 


1524 


.3^2 


2.199 


6677 


800 


7620 


.3^2 


2.161 


13354 


700 


15240 


.3^2 


2.215 


20032 


800 


22860 


.421 


2.705 


4109 


523 


15240 


.421 


2.708 


6163 


521 


22860 


* Based 
gm/cm. 


upon a : 
, sec 


fixed laminar 


viscosity of 


1.78 X 10 "^ 



** Niter is the number of iterations 
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rate, non-dimensional temperature, axial velocity, radial 
velocity, effective viscosity, density, and turbulence 
length scale. Variations of shear stress and heat transfer 
coefficient along the walls were also calculated. The 
location of the flow reattachment point was determined by 
solving for the point at which the wall shear stress was 
zero. 

C. PISTEP II WITH WALL MASS ADDITION 

The boundary conditions were altered along the "north” 
wall as shown in Figure 10. A linearly varying stream 
function (l.e., a constant radial velocity at the wall) 
was used to simulate injection of mass uniformly through the 
wall. For this initial study, the chemical composition of 
the mass injected through the wall was chosen to be identical 
to that of the main stream coming in the inlet. FAR is a 
term that is analogous to a fuel air ratio, and represents 
the fraction of the total incoming inlet stream that is 
injected through the wall. For example, FAR = 0 gives no 
wall mass injection. FAR = 0.5 gives an amount of wall mass 
injection equal to half of the inlet stream so that the amount 
leaving at the exit is 1.5 times the amount entering at 
the inlet. 

FAR was varied in increments of 0.1 for the case in 
which inlet velocity was 500 ft/sec (VIN = 15240 cm/sec) 
for an h/D of 0.342. Stream function output was punched 
onto cards which were used as input for a separate program 
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which generated contours of streamlines using the Calcomp 
plotter (see Appendix A). Figures 13 through 21 show the 
streamline patterns obtained for various conditions. It 
should be noted that the streamline patterns are shown only 
for the upstream half of the ramjet and that the radial 
coordinate has been expanded more than the axial coordinate 
so that the figures are not to scale. This permits a better 
visual presentation of the streamlines in the recirculation 
zone. Note that with wall mass addition (Figures l4 through 
l8), a streamline that originates on the north wall in the 
recirculation region escapes from the recirculation zone and 
enters the main stream of flow. Figures 20, 13, and 19 show 
the effect of increasing step size (h/D) while maintaining 
inlet velocity constant (mass flow rate decreases). Compari- 
son of Figures 20 and 21 shows the effect of decreasing 
inlet velocity (or inlet mass flow) with a fixed h/D and 
comparison of Figures 13 and 21 shows the effect of decreasing 
step size with a fixed inlet mass flow rate. 

Table V summarizes the results for wall mass addition. 

The reattachment point location, based on the point where 
wall shear stress is zero, moves upstream slightly as wall 
blowing is increased. The ’’wall shear stress” is calculated 
from the stream function difference between the wall node 
and the first node out from the wall. In the tube between 
these two nodes the axial velocity is at least an order of 
magnitude greater than the radial velocity. This permitted 
the calculation of the ’’reattachment point” in the same 
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manner as for no blowing. Also, the magnitude of the radial 
velocity at the wall is quite small compared to the magnitude 
of the inlet velocity. It should be emphasized that these 
results were obtained using the same boundary conditions on 
w/r, T, e, and k as for no blowing. Further studies are 
required to determine whether or not this assumption is 
realistic . 



TABLE V 

PISTEP II WITH WALL MASS ADDITION (h/D = 0.3^2) 



Vin (cm/sec) 


FAR 


X /D 
r 




152^10 


0.0 


2.211 


0.0 


152^]0 


0.1 


2.191 


4.7 


152^10 


0.2 


2.168 


9.5 


15240 


0.3 


2.145 


14.2 


15240 


0.4 


2.123 


19.0 


'15240 


0.5 


2.101 


23.7 



Note: 900 iterations for each of the above runs 



D. PMM/AIR RAMJET 

Seventy-four hot firing tests were conducted. A Fortran 
language computer program was written to reduce the data. 

A computer listing together with the output is included in 
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Appendix B. The first three pages of the output tabulate 
the measured data for each test In which combustion was 
sustained. This Is followed by the calculated results. 

This program was modified to reduce the data for the tests 
In which combustion was not sustained. A computer listing 
and output for the non-sustalnlng tests Is Included In 
Appendix C. 

The expected uncertainties In the experimental results 
was calculated using the error analysis method of Kline and 
McCllntock [Ref. 17]. Using the notation defined In the 
ramjet data reduction computer program, the uncertainties 
are given In Table VI. 

The step Inlet caused non-uniform regression of the fuel 
grains. There was less regression of the fuel at the head 
end In the recirculation region, and nearly uniform regres- 
sion downstream of the reattachment point, where regression 
rate was observed to decrease slightly with Increasing 
distance from the head end. The Inlet effect was more 
pronounced for larger values of h/D. For the majority of 
the tests there was a significant amount of carbon deposited 
on the surface of the fuel at the head end. The reattachment 
point locations were measured to the edge of the carbon 
deposits or to the location of maximum regression. These 
two locations were nearly coincident In most cases. The 
regression patterns for three tests are shown In Figure 22. 

Regression rates were calculated using two methods. 

Based upon the exit dlamter, average regression rate can 
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TABLE VI 



ERROR ANALYSIS 





Symbol 


% Error 


Variable 




WDOT 


6.3 


air weight flow rate 




TC 


13.8* 


calculated chamber 
temperature 




RDE 


l.H 


• 

r based on exit diameter 




RWT 


1.1 


r* based on weight loss 




WF 


.6 


fuel weight flow rate 




FG 


8.6* 


gross thrust 




FISP 


8.6* 


fuel specific impulse 




RESTEP 


5.8 


step Reynolds number 




GAIR 


5.7 


air flux rate 




UPORT 


15.0* 


average port velocity 




SFC 


8.6* 


specific fuel consumption 




EQN 


11.0 


• 

empirical r equation 


* The uncertainty in 
the gas constant. 
Complete expansion 


the ratio 
R, were not 
through an 


of specific heats, y> and 
included in this analysis, 
ideal nozzle was assumed 



in calculating FG, FISP, and SFC. 
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be defined by: 



r = - 2 ~^ — (In/sec) 




where d^ and d^ are the initial and final exit diameters 
before and after firing, and At is the burn time. An 
average regression rate based on weight loss of the fuel 
can also be calculated: 



where and are the Initial and final weights of the 

grain, p is the density of the fuel, and L is the grain 
length. This expression for the regression rate gives a 
better average value when there is non-uniform regression 
along the length of the grain, and was used throughout 
this investigation. 

The dependence of the regression rate of the fuel on 

chamber pressure, P (psla) , inlet air temperature, T (°R), 

2 

and average mass flux of air, G (Ibm/in sec), was determined 

( 

by holding two of these parameters at fixed values and 
varying the third for several tests. It was found that 
the regression rate varied as: 




2 At 



(In/sec ) 




(in/sec ) 
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where C is a proportionality constant. For a step h/D of 

- 1 + 

0.333, C = 2.3 X 10 and for a step h/D of 0.416, 

— 4 

C = 3.4 X 10 . A plot of the regression rate based on 

weight loss, » versus this empirical regression rate 
equation is shown in Figure 23. On the sixth page of the 
ramjet data reduction computer program output in Appendix B, 
the numerical results are tabulated. The expected uncer- 
tainty in the empirical equation was calculated to be 11% . 

Muzzy [Ref. 4] has shown that the regression rate 
equation for a hybrid rocket operating in the kinetically 
dominated region can be expressed as: 



f = C 



where a varies from 0.3 to 0.4 and b varies from -0.1 to 
-0.2, and X is the longitudinal coordinate. 

Thus, the average regression rate of the fuel in a solid 
fuel ramjet is very similar to that for pressure sensitive 
hybrid rocket combustion. This is not surprising since the 
reattachment point normally occurred at less than one-fourth 
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of the grain length. Aft of the reattachment point combus- 
tion should approach that of a kinetically controlled hybrid 
combustor. 

The maximum regression rate of the fuel occurred near the 
location of the flow reattachment point. This is due to 
several factors. More oxidizer reaches the wall from the 
impingement of the main oxidizer stream at the reattachment 
point. In addition, the maximum Nusselt number occurs near 
the reattachment point location. 

Data from 36 of the tests in which combustion v/as not 
sustained are tabulated with the computer program in Appen- 
dix C. For the smallest step h/D used (0.25), combustion 
was sustained only while the ignition system was operating. 

As soon as ignition was terminated, the flame blew off. 

For the middle step (h/D = 0.333) the flame blew off as 
soon as the air flow was activated for many of the tests 
at higher flow rates. Some of the tests in which combus- 
tion was not sustained may have been due to insufficient 
Ignition prior to activating the air flow. 

Tests performed at the United Technology Center [Ref. 8] 
have indicated that inlet configuration is more important 
in determining when combustion will be sustained. A critical 
step h/D ratio was found to be the best criterion. For 
values of h/D less than the critical value, combustion could 
not be sustained. It was reported that for one particular 
all- hydrocarbon fuel the critical h/D was 0.1^1. For 
polymethylmethacrylate they found that the critical h/D 
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was 0.31* This criterion did not vary with step Reynolds 
number. 

The data from this investigation indicated that for 
high inlet velocities with h/D greater than the critical 
value determined by UTC, combustion may not be sustained. 
It appears that not only is h/D important, but inlet 
velocity, ignition time, and inlet temperature may also 

be important. At identical air mass flow rates, some of 

\ 

the tests at higher inlet temperatures failed to sustain 
combustion. This is probably a velocity effect since at 
constant mass flow rates, increasing temperature Increases 
the velocity. 

The Pistep II results, although for non-burning cases, 

f 

should indicate trends to be expected in the burning tests 
The data indicate that the magnitude of the velocity at a 
point in the recirculation zone is directly proportional 
to the inlet velocity. Also the velocity in the recircu- 
lation zone V, was found to be Inversely proportional to 

h/D. Thus V « VIN and V “ l/(h/D). It was also found 
■ ^r h 

that p but did not vary with inlet velocity. 

A simplified model of the flame stabilization region 
is one in which the recirculation zone is considered as an 
ignition zone whose length, L, corresponds to the reattach 
ment point location. At the flame blowoff condition, the 
gases traveling through the recirculation zone have an 
average contact time equivalent to the ignition time of 
the gases: t = t.. If V is a characteristic velocity 

C 1 




m 



w 

. 4 




■M 











in the recirculation zone, then t « L/V. Since 

c 



V cc 



VIN 

h/D 



and 



L ^ X /D cc h/D 
r 

then one possible expression for contact time is 

(h)2 

t c. 
c VIN 

The step h/D has a more dominant effect in determining 
flammability limits, however inlet velocity is also impor- 
tant. This treatment is oversimplified and was discussed 
only to indicate that the inlet velocity may be of signifi- 
cance in determining blowoff limits. In actuality, ignition 
time is a function of the fuel-air ratio in the recirculation 
zone and this is in turn a function of h/D. Thus, the func- 
tional form becomes more complex as one considers the 
details of the recirculation region. 

It should also be noted that the simplified model 
described above neglects the effects of turbulence intensity 
on flame stabilization. Examination of Plstep II output 
indicates that for a constant inlet mass flow rate, increasing 
step height Increases both the turbulent kinetic energy and 



k2 



the turbulence dissipation rate. Also for a fixed inlet 
velocity, increasing step height decreases the recirculation 
zone velocity, the turbulent kinetic energy and the turbulence 
dissipation rate. 

There are three blow-off situations to be considered. 

The first one is a steady-state blowoff limit that occurs 
when the mass flow rate in the recirculation zone is greater 
than the mass consumption rate by chemical reaction. The 
mass flow rate through the zone will depend upon the size 
of the zone and the velocity through the zone. This may be 
important in throttling applications during normal combus- 
tion. The second situation is one in which there are igni- 
tion transient blowoff limits during starting. Low chamber 
pressure during ignition gives higher inlet velocities which 
may cause the flame to blowoff even though combustion would 
sustain if ignition occurred at the expected operating 
pressure. Thirdly, it is possible that the flow at the 
step inlet may become choked at high mass flow rates when 
the chamber pressure is low during ignition. Subsequent 
expansion could modify the size and flow characteristics 
of the recirculation zone and cause the flame to blowoff. 

This indicates that if a ramjet operating with high flow 
rates could be started at a higher chamber pressure, 
combustion may be sustained. 

The Plstep II computer program should prove to be 
advantageous in studying blowoff limits. When chemical 
reaction and wall mass injection are Included in the 
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program, blowoff velocities in the recirculation zone could 
be realistically studied. The effect of h/D, inlet velocity 
(or mass flow rate), and chemical reaction rate could be 
considered. Blowoff limits have been calculated in this 
manner for premixed gases in a dump combustor [Ref. 15]. 

If the results of Abbott and Kline [Ref. 1] are similar 
to what may be expected in an axisymmetric situation with 
combustipn, inlet turbulence would not be expected to affect 
the recirculation zone. However, inlet turbulence may 
affect the regression rate of the fuel downstream of the 
reattachment point. Further tests are needed to clarify 
the affects of inlet turbulence and distortion on solid 
fuel ramjet combustion. 



VII. CONCLUSIONS 



The cold flow visualization data were in general agree- 
ment with previous studies. Reynolds number had no measur- 
able effect on the size of the recirculation zone and the 
reattachment point location. The reattachment point moved 
further downstream from the face of the step as step height 
Increased. 

The Pistep II computer program appears to predict with 
acceptable accuracy the reattachment point locations for 
non-reacting flows with no wall mass addition. Pistep II, 
modified for wall mass addition, indicated that the size of 
the recirculation zone and the reattachment point location 
are not greatly affected by wall mass addition. The location 
of the reattachment point moved slightly closer to the step 
face as wall blowing was increased. Experiments are required 
to verify this result, and the assumptions concerned with 
the boundary conditions used in the computer program require 
further investigation. The Pistep II computer program has 
proven to be accurate enough to warrant further development 
as a ballistics program for the recirculation region. 

Solid fuel ramjets which use PMM as the fuel have an 
average regression rate that follows the expression: 
r = C P T G In addition, the maximum 

regression occurred near the reattachment point location. 
Downstream of the reattachment point the regression rate 
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decreased slightly with Increasing distance from the head 
end. Thus, solid fuel ramjets have an average regression 
rate that closely follows the theoretical expression derived 
for klnetically controlled hybrid rocket combustion. 

The step h/D is the dominant parameter in determining 
flame stability limits; however, inlet velocity is also an 
important parameter. 

The results of this investigation indicate that a 
reasonable internal ballistics model for solid fuel ramjets 
may be obtained by using the Plstep II computer program 
(modified for chemical reactions and wall mass addition) 
together with the regression rate equation for klnetically 
controlled hybrid rocket combustion. The modified Pistep 
II program should be used to analyze the recirculation 
region (blowoff limits, etc.) and to provide the input 
data (chemical composition, etc.) to the regression rate 
equation for klnetically controlled hybrid rocket combus- 
tion. The latter . should be employed downstream of the 
reattachment point. 



VIII. SUGGESTIONS FOR FUTURE WORK 



This work was an initial study of the internal ballis- 
tics of solid fuel ramjets. Additional experimental and 
analytical studies are required in order to a) obtain 
increased combustion efficiency, b) refine the internal 
ballistics model and c) study the effect of other inlet ge- 
ometries' and flow conditions on the regression rate of the 
fuel. Specific experimental studies needed include deter- 
mination of a) the effect of wall mass addition on the 
recirculation zone for a step inlet using inert flows (to 
compare with the theoretical predictions using Plstep II), 
b) the effect of different inlet configurations such as 
swirl vanes and plug inlets on fuel regression rate and 
flame stability limits, and c) the effect of different 
mixing techniques on combustion efficiency. Analytical 
studies required are a) modifications to the Pistep II 
computer program to include chemical reaction with heat and 
mass transfer at the wall and b) combination of the recir- 
culation zone model (as input) with the kinetlcally controlled 
hybrid combustion model to effect a usable internal 
ballistics model for solid fuel ramjets. 
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FIGURE 1. COMBUSTION MECHANISMS 
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FIGURE 3. SCHEMATIC OF SOLID FUEL RAMJET APPARATUS 
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FIGURE k. SOLID FUEL RAMJET MOTOR 
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FIGURE 5. RAMJET MOTOR CM TEST STAND 
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FIGURE 6. TEST STAND 






FIGURE TEST STAIJD 





NORTH WALL 



STEP WALL 



^ = 'i'wall " constant 






w/r, k, e 



il) = U) = constant 
wall 

T = 1 

= 0 

r z 

(o/r, k, e calculated from 
wall functions 




V„ = V. = constant 
z in 



ij; from 



t - t 



T = 



inlet 



^wall ^inlet 



(o/r = 0 

k = 0.004 v: 



in 



h / 

c = k^^/l 



I = 0.05 R2(0.22) 
T = 0 



FIGURE 9. BOUNDARY CONDITIONS FOR PISTEP II 
(Figure 2 from Ref. 15) 
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FIGURE 10. BOUNDARY CONDITIONS WITH 
WALL MASS INJECTION 
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COLD PLOW: X^/d = 7.15(h/D) +0.35 



PISTEP 



X^/d = 5.5^(h/D) + 0.36 



FIGURE 11. REATTACHMENT POINT VS. STEP HEIGHT 
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h/D 

K — K Krall and Sparrow - based on location of maximum 
heat transfer coefficient (axlsymmetric flow) 

A — A Abbott and Kline - note different reattachment point 
locations for each side of two-dimensional channel 

C — C COLD FLOW RESULTS 

P — P PISTEP II RESULTS 
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FIGURE 13. STREAMLINES FOR NON-REACTING FLOW 
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FIGURE 14. STREAMLINES FOR NON-REACTING FLOW 
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FIGURE 15. STREAMLINES FOR NON-REACTING FLOW 
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FIGURE 16. STREAMLINES FOR NON-REACTING FLOW 
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FIGURE 17. STREAMLINES FOR NON-REACTING FLOW 
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FIGURE 18. STREAMLINES FOR NON-REACTING FLOW 
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FIGURE 19. STREAMLINES FOR NON-REACTING FLOW 
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FIGURE 20. STREAMLINES FOR NON-REACTING FLOW 
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FIGURE 21. STREAMLINES FOR NON-REACTING FLOW 
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FIGURE 22. REGRESSION PATTERNS 
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FIGURE 23 . EMPIRICAL REGRESSION RATE EQUATION 
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SOLID FUEL RAMJET DATA REDUCTION, NON-SUSTAINING TESTS 
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CALCULATED RESULTS 
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